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ABSTRACT 

The  hydrodynamic  problem  of  the  rotating  disc  electrode  la 
approximated  using  the  weighted  residual  technique  of  orthogonal 
collocation.  Results  are  given  for  a  potential  step  applied  to 
the  electrode  where  simple  reversible  charge  transfer  to  a  single 
species  takes  place. 


ORTHOGONAL  COLLOCATION  SIMULATION  OF  THE 
ROTATING  DISK  ELECTRODE 
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Analytical  solutions  for  ths  current  at  rotating  disc 
electridts  ( RDE)  (Figure  1)  are  known  Cor  only  the  simplest 
electrochemical  reaction  mechanisms  or  electrical 
perturbations.  Orthogonal  collocation  solution  techniques,  like 
finite  difference  approaches,  give  a  single  mathematical  net hod 
tor  approximating  virtually  all  electrochemical  responses  and 
mechanises,  the  ROE  is  described  by  a  partial  differential 
equation  containing  square  distance  terms,  it  is  shown  that  such 
non-linearity  presents  no  special  problems  in  the  approximation 
technique.  The  results  for  ths  ROE  are  available  from  changing 
only  2  program  staps  in  the  general  program  solving  a  simple 
electron  transfer  mechanic  at  a  stationary  planar  electrode. 

DIMENSIONAL  ANALYSIS 

The  general  differential  equation  describing  mass  transport 
to  a  rotating  disc  electrode  is  given  by 


4ci 

v‘  ,f 


where  the  are  the  concentrations  of  electroactive  species,  t 
is  the  time,  g  is  the  distance  coordinate  normal  to  the  electrode 
surface,  Oj  are  the  diffusion  coefficient,  r  is  the  distance 
coordinate  parallel  to  the  electrode  surface,  the  are  velocity 


the  z  axis. 

ft 

In  the  assumption  ot  an  ideal  hydr^namic  model,  there  is  no 
change  in  concentration  in  the  •  direction.  Also,  we  assume  that 
Vx  is  independent  of  r,  which  implies  a  uniform  flow,  tor  can 
thus  approximate  ic/ir  ■  0  over  the  whole  electrode  when  the 
radius  of  the  electrode  is  small  with  respact  to  the  sheath. 

The  assumptions  reduce  the  differential  equation  to  the 

form: 
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The  velocity  term  Vx  has  been  given  by  Uvitch  (1)  and  is  a 
function  of  the  kinematic  viscosity  v  and  the  rotation  velocity 


-0.51  **/2v-l/2#2  .  Kx2  13J 
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Me  require  to  have  the  equation  in  dimensionless  form*  to 
choose  as  the  distance  coordinate 
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The  values  of  each  of  these  matrix  elements  ere  easily 
calculated  by  Mtrix  Multiplication  and  inverse  matrix  routines 
(such  as  those  based  on  Causs-Jordan  elimination)  (41. 

Substitution  of  the  discretized  differentials  into  111) 
yields; 


« 

I  B4.  c  (2  ,T) 
)-l  °  3 


,  N*2 

•V  1  Ail 
*  j-1  13 


(18) 


We  have  now  M*2  equations  (Z«0  and  2-1  plus  the  N  central  roots) 
in  N*2  unknown  c'iZj.T).  The  dinension  is  reduced  by  two,  and 
the  problem  specified  to  final  fore  by  introducing  into  (18)  the 
known  boundary  conditions 
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which  are  *he  first  and  laat  terns  of  each  sunnationi 
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dc 
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to  yield 
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r  B4,c  (2  ,T)  - 

j-2  13  3 


)A 

Thus  the  concentration  profiles  as  a  function  of  tine  are 
readily  and  efficiently  generated  by  Integrating  the  N  sets  of 
first  order  differential  equation  |21|«  One  may  use  available 
subroutines  fro*  software  libraries  such  as  IHSL  or  SSP,  or  nora 
specialized  programs  taking  the  stiffness  of  the  differential 
equations  into  account  (2 >4). 

As  has  been  shown  previously#  More  complicated  mechanisms 
My  be  readily  Incorporated  into  the  equations  (3,4)  *•  well  as  a 
variety  of  electrochemical  techniques,  electrode  geometries 
( S-8) ,  and  edge  effect  considerations  (9,10). 
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ROE  CURRENT 


Current  to  the  roe  is  given  by  the  flux  at  the  surface i 
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or  in  terns  of  the  distension  leas  concentration  [81  and  distance 
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which  is  given  immediately  during  the  simulation  at  each  tine 
point  by 
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since  the  tan  in  brackets  is  the  diaens ionless  flux  in  tqt  ion 
(24|  at  Z«0  by  equation  «ii|. 

Ne  can  compare  our  result  to  a  series  solution  ()  which 
states  that  the  time  dependent  current  is  proportional  to 
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Introducing  the  dimensionless  time  variable  M  intofcM  we 
find  the  current  is  proportional  to 

1  ♦  *  £  expl-jW  R2/2T)  1201 
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,  so  that  the  ratio  of  the  current 
calculated  from  equation  (251  to  the  expression  1 281  should  be 
constant.  Results  for  s  given  set  of  parameters  are  given  in 
Table  1.  it  le  seen  that  for  only  ala  collocation  points,  the 
current  values  are  eccurate  to  0.05%*  This  value  aay  be  further 
decreased  by  forcing  higher  accuracy  in  the  integrations  (2)  or 
increasing  the  number  of  collocation  points,  although  ten  points 
do  not  exhibit  significant  improvement  over  nine.  The  maximum 
relative  error  decrease  using  ten  instead  of  nine  collocation 
points  is  0*0071. 

The  results  for  the  rotating  ring  disc  electrode  are  being 
optimised,  and  will  be  the  subject  of  a  forthcoming  publication. 
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T  Simulated  Flux 
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10'* 

0.7*6 

10"4 

0.7702 

10' 3 

0.7700 

10' 2 

0.7700 

10" 1 

0.7700 

10° 

0.7704 

Table  1.  Ratio  of  simulated  flux  to  analytical 
current- time  function  (Ref.  11)  Bata  •  10  1. 
For  typical  rotation  speeds  (1000  RPM), 
viscosities,  and  diffusion  coefficients,  the 
real  total  time  simulated  in  the  table  is 
t  s  100  ms. 
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